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INTRODUCTION 


This project applied modern theories of turbulence and fluid mechanics 
to astrophysical phenomena. Of special interest were periodic events on 
the surface of the sun. Research and papers written during the grant are 
summarized below. 

1 . stellar structures 

Streamers in the solar corona are thought to be caused by structures 
in the underlying magnetic field. However, we have found that most expanding 
gases will form streamers and expand non-uniforroly , There roust be a retarding 
force (such as gravity in the case of the solar wind) and the expansion must 
be subsonic. 

Our work on this problem was presented at an APS meeting. The abstract 
is enclosed as Appendix A, The paper titled "Natural Inhoroogeneities in an 
Expanding Gas" is Appendix B. We showed that the expanding gas had a linear 
instability to initiate streamer formation, but were unable to use turbulence 
to find a non-linear limiting process for the instability. 


2, TURBULENCE IN A MAGNETIZED PLASMA 

The Orr-Sommerfeld equation governs growing instabilities in a fluid 
boundary layer. Its solutions determine the Reynolds number at which 
turbulence develops. We have added a magnetic field and electrical resistivity. 



This problem applies tc plasma physics (fusion research) and to several 
astrophysical situations. In particular, the interior of the sun has 
boundary layers between regions of magnetic plasma flowing in opposite 
directions. It Is possible that turbulence In these boundary layers 
could have a characteristic frequency that leads to observable phenomena 
on the sun’s surface. 

Appendix C is titled "Onset of Turbulence in a Magnetized Boundary Layer." 
We found that a higher Reynolds number is needed to start turbulence in the 
presence of a magnetic ■‘■"’eld because energy is required to bend the field 
lines attached to the fluid. 


3. cosmological TURBULENCE 

Appendix D is titled "Phase Coherent Turbulence in the Early Cosmos." 
Turbulence in fluids has prefered frequencies and eddy sizes. If the 
cosmological gas was turbulent shortly after the big bang, then galaxies could 
have been formed by turbulent eddies. This hypothesis is still open because 
the Reynolds number is very high in astrophysics. It is not clear how to 
extrapolate from laboratory experiments at Re ■< 10 to astrophysics at 
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Re >10^^. 
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Abstract Submitted 


for the Washington fleeting of the 


American Physical Society 
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Physical Review Bulletin Subject Heading 

Analytic Subject Index in which Paper should be placed 

Number — 95j| Astrophysics 

Bulletin of the American Physical Society 597 (1980). 


Coronal Structures without a Magnetic Pield. * 

M. A. CROSS, Gr ambling U . and J. A. JOHNSON III, 

Rutgers U . — It is believed that structures in the 
atmosphere of the sun are controlled by the background 
magnatic field. We show that coronal streamers can be 
formed without a magnetic field. If theta dependence 
is kept in the Navier-Stokes equations for the solar 
wind, then a density enhancement will grow. This growth 
is followed in the non'"linear equations until a streamer 
is formed. 'Viscosity stops the streamer's growth when 
there is a large difference in speeds inside and outside 
of the streamer. The viscosity needed to reproduce 
observed streamers is compared with the classical 
viscosity. 

*Supported in part by NASA grant NAG-8002 


Submitted by 
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Rutgers University 
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Joseph A. Johnson III 
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ABSTRACT 

Using classical fluid mechanics and a latitude dependent 
hydrodynamical model, we find that unmagnetized perturbed flow 
evolves into high and low density regions. The growth 
mechanisms for density enhancements are discussed along with a 
new nonlinear numerical solution for their large amplitude 
development. 

*Work supported in part by NASA grant NAG 8002 and NSF grant 


SER-810430. 



Examples of expanding gases occur in many astrophysical 
events, e.g, the Big Bang, novae, stellar winds, and planetary 
atmospheres. These processes are usually modeled with a 
restriction to spherically symmetric hydrodynamics. It is 
impossible for such models to produce structures whose 
parameters depend on latitude and longitude. The structures 
actually formed in these expanding gases range from galaxies on 
the cosmic scale down to streamers and bursts on smaller 
scales. Some of these structures might be formed by 
hydrodynamic processes alone. In this paper we study the 
characteristic structures that might evolve from the 
unrestricted Navier-Stokes equations. Two hypotheses are 
investigated: "The expanding gas would still have structures 

even if there were no magnetic field," and "Any subsonic 
expanding gas might break up into structures." 

Stellar streamers are now thought to arise from the regions 
in stellar atmosphere where collisional phenomena predominate 
(Kaplan, Pikelner, and Tystovich, 1974). In such regions, a 
hydrodynamical approach is appropriate for the analysis of the 
evolution of density perturbations (Lemaire 1978). This 
approach should be particularly relevant to those stars with 
weak fields or to stellar atmospheres for which the onset -of 
amplified fluctuations is only weaky dependent on the strength 
of the local magnetic field. By way of illustration, the 
calculation will be applied specifically to the solar corona 
because its parameters are well known and because it has a 
subsonic region. We shall use . solar parameters, even though 
the model is not applicable to the solar corona 
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where bhe magnetic field dominates. Of course/ the parameters 


c,m easily be changed as new data become avilable to apply to 

OrtlGlIMAL 
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The existence of streamers in an expanding gas shows that 


other astrophysical situations. 


spherical symmetry is broken when random density perturbations 
grow and become stable configurations. Ignoring explicit time 
dependence (as is the conventicn')/ the latitude dependent 
Navier-Stokes equations are 


CnvsInB) = D V' 


n 


( 1 ) 


CvoK'T) - ^ (V- ^ ^ Vjt) 


ni. _ _L + (VI) 

~ r ' Yr\r d-Q’ 


( 2 ) 


(3 


±.kuli . - is- -rM), (rc)^vi) 


(4) 


where: n is the total particle density; m is the average mass 
per particle; u is the radial speed; v is the latitude 
component of velocity ; (TC) and (yi') symbolize thermal conductivity 
and viscosity terms respectively; and the transport 

coefficients are (Scarf and Noble, 1965) yK = 8.3 x 

-17 , , -9„^373, 

10 '■y gm/cm-sec, D = h/mn, K = 6 x 10. 

erg/cm-sec-K for qur purposes. Details of the thermal 
conductivity and viscosity terms may be found in Yuan (1967) < 
Linear Treatment . (See Cross and Blockwood 1974.) The zero 
order unperturbed solution ignores theta dependent terms in 
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eqs. (X) - (4). We add a perturbation (indicated by the 
subscript 1) which peaks at angle 8*^ with width “Z : at 

= - (nj /X^), 

( "V n/&^) = 0 with similar forma for T and u perturbations. 

The fourth function v is zero at 8*^ but has nonzero 
since the perturbed pressure (from n^ and causes 
sideways expansion away from 0-^. For a rapidly growing 
perturbation, we assume a small radial scale length, i.e. 
(,>nj^/lfr) = o^n^i^ yy (n^^/r). Equations (1) - (4) in linear 

form without viscosity, diffusion, or thermal conductivity are: 

+ o< n Cn/r) ^ 

YY\ Yi u.ai Vlj A' ~ }< ot ( Y)^T i n7^ ) 

n U C8v/^(&) C K/m r ( 51, T'ivjT’ ) 
n, i - ( VI 77 

respectively. The gravitational term which would have appeared 
in eq. (6, has been dropped because it is much smaller than the 
pressure term; thermal energy is comparable to gravitational 
energy (kT js:-GMir,/r) but the perturbed pressure should have a 
steep gradient. -That is, (dn^^/dr) = 0 ^ n^ >■)> (n^^/r). 

Eliminating T^ and both components of velocity from eqs, (5) 

- (8), and assuming mu <4,5kT/3, produces c^= 1/rt: or , 

n,(r,9-) = n^ (0*) , Since is the width of the 
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( 6 ) 

(7) 


perturbation, narrow perturbations or streamers will grown 
faster as they are convected outward by a stellar wind. 

Physical Mechanism . The density perturbation grows because 
of interactions between eqs, (5) - (7). 'Ignoring the 
temperature change in an isothermal corona reduces the radial 
momentum eq. (6) to 

, Jil l^T I dn, ORIGSNAiL PACS ES' 

"IZ ^7? OF POOR QUALHY (9) 

that. is, fractional change in velocity is proportional to the 
fractional change in density. Recall that kT mu below 
about 5 Rg ; thus velocity changes have the larger magnitude 
below 5 Rq* Eq. (9) also shows that changes in 
density make a pressure gradient which drives the velocity. 

When the continuity equation is applied to a small volume 

at a fixed position in the corona, one sees a density 

enhancement convected upward into the volume by the stellar 

wind. This makes a locally higher pressure which is smoothed 

out by sidev/ays expansion out of the volume governed by eq, {'i) . 

As the density enhancement passes through the volume, it loses 

2 

flux due to sideways expansion. The particle flux nur 
becomes successively smaller after passing through volume 
elements at successively larger r. The only way to both 
decrease the flux and also satisfy eq. (9) is for the speed to 
decrease and the density to increase, i.e., for the density 
enhancement to grown. 
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Nonlinear Numerical Solufcion, Models of expanding gases 


usually assume initial conditions at one radius and then eqs. 

(2), and (4) are numerically integrated to solve for n, u 
and T at other radii. We add eg, (3) and the 0- component of 
velocity, V, to this process. 

Recall that the solar wind has a critical point at the 

supersonic transition where 5KT/3 m u. . Whang, Liu and 
» 

Chang (1966) have shown that the singularity can* be avoided by 
retaining the second order terms in the viscosity terms 

of eg, (2) and in the thermal conductivity terms in 

eg. (4). For this reason the numerical work reported here was 
done with eqs. (2) and (4) rearranged to use second order 
radial derivatives. 

Equation (1) was left at first order in by ignoring 

diffusion on the right side. Finally, in eq, (3) it was 
assumed that the viscous term was unimportant and 

could be approximated numerically from changes in at 

previous (lower) radii. Therefore eq. (3) was kept in the 
first order form for 9-V/^r as written above. 

The functions of were calculated on grids of 30 to 180 
points equally spaced in latitude from 0° to 90°. Symmetry 
about the equator (at 90°) was assumed. The original 
pressure was preserved by using both a density perturbation 
n^and a temperature perturbation T^/t^ = -n^/n^. The 
numerical solutions used eqs. . (i) -* (4) to calculate the 
radial derivatives at all points in the theta grid at a 

- 6 - 


OF POOR QUAinT 


OF POOR QUALITY 

starting radius. Then these derivatives were used to find the 
functions at a higher radius. V/e made the usual check of 
changing the radial step size; starting from a large radial 
stePf the step was decreased until the solution was independent 
of step size. 

The perturbation was a Gaussian in latitude- Other initial 
conditions were typical of solar wind models. Figure 1 shows 
the total praticle density plotted three times as the solutions 
precede outward in radius. We see that a density perturbation 
grows as it is convected outward. Figure 2 shows the history 
of the amplitude of the central maximum relative to the two 
minima on either side. Data from the numerical solution are 
plotted in this figure and compared favorably with the curve 
(amplitude)^ P • with X. - 5.7°. Notice that the 5.7° 
wide Gaussian would also be a reasonable fit to the peak in 
Figure 1 (c). Since the initial perturbation has a Gaussian 
form^ the short wavelengths in a Fourier decomposition of the 
perturbation should predominate and give the growing structure 
an increasingly narrower width, (See the linear growth rate 
derived above.) The results shown in Figures 1 and 3 confirm 
the expected behavior with the peaks getting progressively 
narrower and the distance between the centrc.1 maximum and 'the 
two minima on either side gradually decreasing. 

VJe have not determined the ultimate nonlinear development 
of these radial streamers. Neither have we provided a 
mechanism which* prevents their angular widths from decreasing 
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indefinitely until they are unobservable* However, it is easy 
to speculate that a mfjfjnetic field and turbulence would 
complete the model which has just been introduced. The main 
effect of an imbedded magnetic field on the instability 
presented here would be that of making it difficult for 
particles to cross field lines. This would limit the 
transverse velocity and slow the growth of the instability. 
Furthermore, this convective instability produces neighboring 
streamers with very different radial speeds. The high rate of 
shear should lead to a transition to turbulence. The effective 
transport coefficients would then increase above the laminar 
flow values and thereby limit the ultimate evolution of the 
amplitude and width of a streamer. 

Nonetheless, we have shown that the perturbed expanding gas 
might evolve into altejrnating streams of high density/low speed 
and low density/high speed regions from purely hydrodynamical 
considerations. This onset mechanism would apply to any 
appropriate ast rophysical environment to the extent that 
hydrodynamic boundary conditions are . applicable in an 
underlying collisional atmosphere. Wherever the expanding gas 
is subsonic, perturbations will grow into radially oriented 
structures. 
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Figure Captions 


Figure 1. The evolution of the latitude profile of a 
streamer-like densiuy perturbation, (a) starting 
conditions: At r « 1.3 x 10^^ cm/ T = 1,45 x 10^ K, u = 

1,4 X 10^ cm/sec/ n = 5,1 x 10*^ cm"^/ dT/dr « -2 x 
10“^ K/cm/''and A n « (0,l)n exp£{ V(/o'’ f ] (b) 

Density profile at r = 1.4 x 10^^ cm, (c) Density 
profile at ,r - 1,66 x 10^^ cm. The particle density is 
n, the solar latitude is fr and r is radial distance drom 
the center of the sun. The smooth curves are the results 
from the numerical solution. 

Figure 2. The evolution of the amplitude of a solar 

streamer. 4n/n =, (n„^„-n„,. )/(0. 5) (n„^„+n„. ) . 

' ' max min ' ' max min 

The results from the numerical solution are given by the 
symbol x; the dashed curve is a plot of A n/n = (0. 00 , 22 ) 

( r/lO^^Jf where r = 0.1 rad («5.7°). 

Figure 3. The evolution of the width of a streamer. The 

width/ / is defined as the angular width at half of th 

central maximum in the perturbation. 
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ABSTRACT 

When the Orr-Sommerf eld equation is derived for a conducting 
fluid with an imbedded magnetic. field , the first appearance of 
amplified perturbations occurs at a higher Reynolds number than is 
expected without the field. 

*Work Supported_ in part by NASA Grant Nag 8002 and by NSF Grant 


SER-810430. 




0R!C21WAI, P’A©?. 

, OF POOR QUALITY 

Solutions for th**- Or r~Soinmerf eld equation which determines the 

stability of small disturbances in neutral laminar flow past a wall 

are well known^. Above a critical Reynolds number the 

disturbances grow and eventually develop into turbulence. However, 

solutions for the analogous problem of boundary layer flow in a 

magnetic field have not been achieved, in this case there should be 

t 

important implications for both astrophysical^ and fusion^ • 

( 

research since, at the present time, the impact of turbulence on the 
transport properties of magnetized media is not well understood. In 
this paper we have extended the Orr-Sommerf eld equation to the case 
of a weakly perturbed magnetic field in a resistive Blasius boundary 
layer. 

The coordinate system is shown in Figure 1. The field is 
parallel to the flow so it has no effect on the unperturbed flow. 

We derive the new Or r~Sommerf eld equation using standard notation^ 
where the zeroth and first order functions and the stream function 
are defined by 


=(u^(y), 

(x,y) = 


0, 0); Bq= (Bp, 0, 0 ); 

0); tP (Y)£ 


bi()t/y) = (bj^ 

; «- 


. by- 0) 

- lii’ .v- 


( 1 ) 


in which pi, and are the wavenumber and frequency of the disturbance. 
If both and u^^ are of the form f(y) exp i( oc x - *^t) then 
the magnetic field equation leads to 

C/?'v)v* - C V 

t by ^y ^0 

C- U.^ J. ol by (2) 
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in which •/(rj«; tr and ^ are respectively the conductivity and 
magnetic permeability. Introducing the stream function, eq. (3) 
becomes 


- .r.7 


X. ^ Py 


^ :=■ by ty -t 


(4) 


From conservation of momentum one finds 


CV3* V -V)^ =• +pV^V 


(5) 


Taking the first order x and y components of (5). 




r 


/ 




P, -f (^X ^'7 " \ 4x ) "^y ^ ^ 


By eliminating the pressure between (6) and (7) and using the stream 
function, one achieves 


if +A c<^ T t <xqj 

r? 2 / V 2 *^ 

_ _ i!l (_0(. ty Px) (-rt > + -V J 

f"- 


(8) 
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The dimensionless variables are; | 

^ Id . ^ /J ^ /d } (^o i 

p, ^ U^J/u ; ^ -=^ ^ Av„cl ; - Uc/u^ 

A/^ - 

‘ » 

where d is the width of the boundary layer, is the magnetic 

Reynolds number, N. is the Alfven number, and U is the fluid 

^ m 

velocity at y =o<a . We drop the use of bars and assume that all 
symbols after this refer to dimensionless variables. Equations (4) 
and (8) become 


^3 h ^ ^ ,) xlci jo -b rk 

I 1 ' t / ,"» y - 




(10) 


"'X C ^ LUo"^ ^ 4o< ^ (iia) 

- ^ / /V/) ^ -U.C. t y by /^ ) 

(lib) !; 

where the prime notation means ^ , 

The left side of eq. (11) is the usual Orr-Sommerf eld equation. 

The right side is our new result. It shows that the first order 
magnetic perturbations in the boundary layer approximation simply I 

add a new term' that depends on the Alfven number and the magnetic i 

Reynolds number. I 
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As Rj^-^0 there is no new effect because the low conductivity 
permits pr»rticles to freely flow across magnetic field lines, 

in order to solve these equations, we have altered the method of 
Mack^ to include our equations (10) and (11), At y «o^: 


-c-y -cy 

a„too)=-l j 


( 12 ) 


TWO analytical solutions are known in the free stream far from the 
wall. The first, called the inviscid solution, has c= o(, in eq. 
(12). This requires i (ot-^ ) in eq. (10). The second 

"viscous" solution has 4^c><exp (-py) with RCoC—^) 

when there is no field. We find f> when Bq^^O by first eliminating b 

between equations (10). and (11a). The result is linear in iP and 

* 

can be solved to yield a cubic equation for p as a function of 
di r R/ Rj^f and from which p is obtained numerically. 

The boundary conditions at y = 0 are: 


Mack's method of solution consists of beginning with an initial 

estimate of oC and -p and then numerically integrating the two known 

solutions from y =o0 toward y = 0. The eigenvalues ot, and are 

altered until a linear combination of the two solutions fits the y = 

0 conditions. Our solutions reduce to Mack's in the limits R 0 

m 

and B^O. 
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We have solved for the wave number o< as a function ofy?^R, 

Rj^, and N^. (This is called the space-amplified case, i.e, a 
perturbation at frequency 7^ grows spatially at the rate exp (ic<x) 
as it is convected along the boundary.) The dimensionless frequency 
parameter F - f /R remains constant as a perturbation is convected 
along the boundary layer. (Schlicting^ ) . Figure 2 presents the 
region of instability in the F-R plane for comparison purposes, 
inside the curves, o<^ is negative and perturbations will grow. 

X'he curve for Rj^^ = 0 reduces to the non-magnetic case investigated 
by others. 

Next we consider changes in the magnetic energy density measured 

by N^. When b^<<Cf , eg. (lib) is linear in tp and the 

quantity reduces to a single parameter. In Table I 

we keep this parameter constant. The eigenvalue is nearly 

constant for the entire range. The growth rate c><- is constant for 

In Figure 3 we compare a set of critical Reynolds numbers 

(that is, the minimum Reynolds number at which perturbations will be 

amplified) with a corresponding set of magnetic Reynolds numbers. 

Notice that the curves of Figure 2 shift toward higher R with 

increasing conductivity and the curves of Figure 3 show that the 

flow will require a new and larger R before perturbations are 

amplified. For example: given the fluid parameters R^^ = 1 and R 

= 1300, turbulence could not develop if the magnetic energy density 

2 2 

were greater than B //.;v .1. 

Thus, as the' energy in the magnetic field approaches the kinetic 
energy of the flow, its influence becomes increasingly important 
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through the inhibiting effect of the magnetic field on the 
$ 

amplication of stochastic fluctuations. Specifically, in the 
presence of a magnetic field, a conducting fluid experiences an 
increasing delay (i.e., the critical Reynolds number increases) i 
the onset of turbulence. 
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FIGURE CAPTIONS 


Figure 1. Coordinate system for the Orr~Sommerfeld Equation. A 
viscous fluid flows past a semi-infinite wall. The 
fluid has an imbedded magnetic field. 


Figure 2. Neutral stability curves for disturbances in magnetic 
laminar flows. F « -^/R is the frequency of the 
disturbance. The imaginary part of the wavenumber is 
negative inside the curves indicating that 
disturbances will grow in space. The changed values of 
the magnetic Reynolds number, R^^^ are: a;aO;b, 


Figure 3, Critical Reynolds number (R ) for turbulent onset as 
a function of the magnetic Reynolds number (R^^) , 

The leftmost curve, a, shows the case for magnetic 
energy density, (B /yt^jpv ), 10 times the kinetic 
energy density; the center curve b is for magnetic 
energy equal to kinetic energy; and for the right 
curve, c, the magnetic energy is one tenth of kinetic 
energy. 


pp. 456-488 (1979). 
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Table Caption 


Table I. Sensitivity of Solution Eigenvalues to Changes in 
Magnetic Energy Density., Herei R « 1280; « .0768; F-Wr « 

6 X 10 and ' 
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TABLE I 



m 

®<r 


1000 

.0001 

.228487 

.177975 X 10 

100 

.001 

.228487 

.1776 X 10 “^ 

10 

.01 

.228487 

.1739 X 10 “^ 

5 

.02 

.228485 

.16976 X 10 “ 

1 

.1 

.22848 

.1368 X 10 "^ 

.5 

.2 

.228474 

.97 X 10 “^ 

.2 

.5 

.228474 

-.3086 X 10 " 

.1 

.1 ■ 

.22842 

-.2460 X 10 “ 

.01 

1 

.227665 

-.0005188 

.005 

20 

.22662 

-.001248 

0 

0 

.23056 

-.0087454 
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In collisional fluids and plasmas, reaction-diffusion 
instabilities arise because of an inhibition of dispersion by 

I 

a temporally competitive non-equilibrium process . A use of 
the Klimentovich formulation has predicted the recently observed 
reaction distortions in unperturbed bursting nonequilibrium flow 
at transition to turbulence . The use of a reaction* diffusion 
mechanism provides a good qualitative explanation of this 

3 

effect as well as the soliton-like turbulent bursts found in 

4 

collisional shock fronts . In this paper, ^e shall explore 
the implications of this mechanism for the early condensation 
of galactic matter. 

First, -lets reconsider the role which thermal fluctuations 
might have played in the formation of galaxies. Classically, 
one begins with a balance between the gravitational instability 
and pressure oscillations for an inhomogeneity of mass M in a 
universe of total initial chaos. The limits on the strengths 

i 

of the initial inhomogeneities are set by the observed degree of 

. 5 

isotropy in the microwave gackground. A critial mass is derived , M_ 

whose evolution in time is summarized in Figure 1, from v;hich 

is derived a minimium undamped mass not in disagreement with 

present observations. These same arguments can shov; that 

ttie preferred galactic shape should be pancake rather than spherical 

and can give a good qualitative relationship between the final 

radius of a gal'axy, its total original mass^ •(•i-n '-the ’beginning) 

and its total angular momentum. 

When- total in,i-tial chaos is removed, the sourcej of instabilities 

0 

in the early cosmos become the same as that of any fluid of compci-t-a Md 
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thermodynamic parameters. Thus, the onset mechanisms are 
determined by the natural evolution of appropriate local 
scales and the development of natural local instabilities 
in the case of the collision dominated epoch, in this context, 
therefore, natural thermal instabilities (rather than the 
artificially imposed hierarchy derived post hoc) can generate 
turbulence. Furthermore, if the front of the evolving cosmos 

can be treated as a tangential discontinuity, then the recent 

, *1 
discoveries concerning the application of Orr-Sommerf eld 

techniques to transition at discontinuities with natural 

instabilities can be applied^. 

Consider turbulent eddies v/hich are locally subsonic in the 

radiation era. Their importance to the growth of turbulence 

is determined by the local Reynolds number. If we ignore magnetic 

effects, the relationship betv/een characteristic turbulent time 
•> 

and critical Reynolds number is shown in Figure 2 for typical 
cosmological values. ‘ By comparing- these results with Figure 1, 
we see that the onset of turbulence is delayed substantially beyond 
the classically derived times; the local Reynolds numbers are 
too lov7. Thus, this new model predicts that galactic condensation 
v/as much later than is usually thought, well into the matter-dominated 
era . 

Next consider turbulent eddies which are locally supersonic 
in the matter-dominated era. Since the turbulent intensity is 
still quite -'small, the onset mechanisms are still determined by 

7 

the techniques ‘from, incompressible turbulent theory . At 
transition, the energy of the turbulence iSj^stitirely contained in 
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the phase-coherent turbulent bursts. Hov/ever, since the universe 

is continuing to expand, nonturbulent regions arehxper i rnclng 

a continuous decrease in local Reynolds nuiubni;' and ai'’; ihorcforc 

increasingly unlikely to produce turbulent f lucbuatinnn v.'hich 

are amplified. This forces a reinterpretation of the meaning 

of a new Mj, viz,,, it should represent a maximum as well as a 

minimum on the typical galactic size, in accordance with observation. 

However, the condensing turbulent bursts interact with the surrounding 

sea of gases; these interactions produce, due to the naturally 

probabilistic characteristics, a variety of inhomogenieties . 

This variety, in turn, leads to the taxonomy of galaxies ranging 

from the elliptical to the various spiral types to the irregular 

shapes. Nonetheless, the frozen in feature of the spiral arms 
c\ 

is^resonable consequence of evolution with minimum surrounding 


cea interaction and fixed initial turbulent parameters. 

~T tv d 

.new consequences can be determined. (1) The class 
or spectrum of galaxies which evolved during the early epoch must 
be complete. New galaxies will, if they are produced , result 
f romjentirely different physical phenomena. (2) If the pre- 
sently observed galaxies all developed during the early epoch 
at or near the critical Reynolds number, then an upper limit 
on the "missing mass" of the universe can be calculated. That 
is, in this model given the size of a typical galaxy, the total 
mass of the universe is no longer a free parameter. ( 
^hese-cnset-mechanisms- imply a distortion of reaction processes, 

ft 

“substantial differences in the distribution of nuclear . matter 

0 

must exist between galatic and intergalatic space. Specific 
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quantitative predictions are being calculated and will bo reported 
in a I'apor to follow, 




